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ABSTRACT 

The late-time optical/radio afterglows of 7-ray bursts (GRBs) are believed to be synchrotron emission of 
electrons accelerated in relativistic collisionless shocks propagating in the ambient medium of the sources. 
However, the fraction / of electrons that are coupled to protons and accelerated remains unclear and a large 
number of thermal electrons that are not coupled to protons may be left behind. If / < 1, the true explosion 
energies of GRBs are f~ times larger than those commonly estimated with / = 1 . Thus the value of / gives 
an important constraint on the nature of the central engine of GRBs and the physics of collisionless shocks. 
Although early-time radio observations can probe the thermal electrons, they are difficult at present. We show 
that the Faraday rotation effects of the thermal electrons may suppress the linear polarization of the afterglow 
at frequencies higher than the absorption frequency in the late time, if the magnetic field is ordered at least in 
parts, and that / can be constrained through the observation of the effects. We find that those effects may be 
detected with late-time, > 1 day, polarimetry with ALMA for a burst occurring within 1 Gpc (i.e., z — 0.2), if 

/- 10 . 

Subject headings: acceleration of particles — gamma rays: bursts — polarization — radio continuum: general 



1. INTRODUCTION 

The afterglows of 7-ray bursts (GRBs) have been observed 
mainly in the optical and radio wavebands in the late time (i.e., 
several hours after the burst trigger) since the late 1990's, and 
they are widely explained as due to synchrotron emission of 
electrons accelerated in relativistic collisionless shocks driven 
into the ambient medium of the GRB sources. The syn- 
chrotron emission mechanism is supported by the detection 
of linear polarization at th e level of ~ 1 — 3% in several opti - 
cal afterglows (for reviews. ICovino et al.l2004HLaz zati 2006). 
These understandings are also being confirmed by recent ob- 
servations with Swift satellite, although the situation is very 
compl icated in the X-ray band and in the early time (for re- 
views, EraillOOl iMUMmillOOl |Zhangl2007|). 

In the standard external shock model of GRB afterglows, 
the late-time dynamics of the shock is determined by the ex- 
plosion energy E and the ambient medium number density 
n. We usually treat the fractions of the explosion energy 
that go into the magnetic field and electrons and the frac- 
tion of electrons that gain the energy of protons as free pa- 
rameters €b, e e , and /, respectively, since these have not been 
derived from the basic principles. We additionally assume 
that all the electrons that gain the proton energy are acceler- 
ated to form a power-law energy distribution, dn/d^ e oc r y~ p 
for j e > 7„, ~ e e (m p /m e )T, where V is the Lorentz factor 
of the shocked fluid, and a thermal component with j e ~ 
e e (m p /m e )T is not produced. We commonly constrain the 
parameters \E, n, e#, e e , p| by the observations, ass uming 
/= 1 (e.g., lPanaitescu & Kuma3l2002tlYost et al.ll2003l) . 

This implies that current observations cannot constrain the 
electro n-proton coupling parameter /. lEichler & Waxmanl 
(2005) showed that the observations also allow the external 
shock model with the parameters chosen as {E 1 = E/f, n' = 
11 1 '/> e B = e B.f, e ' e = e ?/> p' = P] m which the fraction / of 
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total electrons gains the proton energy and the fraction (1 — f) 
is thermal electrons with j e ~ T, as long as m e /m p < f < 1. 
(Hereafter we call the former and latter electron components 
"accelerated electrons" and "thermal electrons", respectively.) 
The parameter / gives an important constraint on the true ex- 
plosion energies of GRBs. It also gives a clue to unveil the 
physics of collisionless shocks. 

The thermal electrons provide additional emission and ab- 
sorption for the afterglow flux only for v < v m , where v m is the 
characteristic synchrotron frequency of the thermal electrons. 
This results in the sharp decline or sharp rise of flux as v m 
passes the observed frequency, which can be detected through 
early-time radio observat ions, t < 10 3 s for v > 10 11 Hz 
(Eichler & Waxman 2005). However, such observations are 
difficult at present. Furthermore, it is unclear w hether the 
standard model is applic able forf < 10 3 s (see e.g. Jloka et al.l 
l2006HToma e"t~aTJl2006l) . 

In this paper, we show that the thermal electrons give the 
Faraday effects on the afterglow polarization substantially 
even in the late time and that / can be constrained through 
the observation of the effects. Those effects may be signifi- 
cant (e ven for v > i> m ) if the magnetic field is order ed to some 
extent (iMatsumiya & Iokall2003l ISagiv et al.ll2004l) . The fre- 
quency below which the Faraday rotation effect is significant 
is expected to be orders of magnitude higher than that sup- 
posed so far with / = 1, because the position angle rotation 
depends on electron Lorentz factor as Ax oc (lnj e )/^. 

2. TRANSFER OF POLARIZED RADIATION 

The transfer of polarized radiation through spatially homo- 
geneous plasma with a weakly anisotropic dielectric tensor 
may b e described by the transfer equation of Stokes param e- 
ters (lSazonovll969l:lJones & O'D^fT^lM^h^fToSOallbh , 
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where s is a length parameter along the ray path, and the right- 
handed system of coordinates with the wavevector k along 
axis 3 and the magnetic field B on plane 2-3 is adopted. Here 
rjj.Q.v are polarization-dependent emissivity, and Kj&y (Kq y) 
are the transfer coefficients related to the anti-Hermitian (Her- 
mitian) part of the dielectric tensor, describing polarization- 
dependent absorption (the polarization of the normal modes 
of the plasma). If \k v \ 3> \k%\, the normal modes are circu- 
larly polarized, and the transfer equation (fl]i indicates that the 
conversion of Q and U occurs. This is well-known Faraday 
rotation. If \k q \ ^> \Ky\, the normal modes are linearly po- 
larized and the conversion of U and V occurs. This is called 
Faraday conversion. 

We define the optical depth r = k/s, the rotation depth 
i~v = \ n v \s, and the conversion depth tq = \kq\s. The proper- 
ties of the solution of the transfer equation ([1} are as follows. 
First, suppose that the absorption effect is not significant, i.e., 
r <C 1 . In this case the equation CP ma y be integrated easily 
dMelrosdll980b1; Uones & 0'Delllll977l) . For \k 
obtain the linear polarization 
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and the circular polarization is given by the intrinsic one, 
lie = \V\/I = \i]v/vi\- F° r T v 3" 1> the linear polarization 
damps. This results from the fact that the emission from dif- 
ferent points through the source have its polarization plane 
rotated at different angles. Analogously, for \k%\ 3> \k v \ and 
tq ^> 1, the damping of lie occurs while ITx remains intrinsic. 

Secondly, in the case in which the absorption effect is sig- 
nificant, i.e., r ^> 1, we may obtain the polarization degrees 
approximately by eliminating the differential term from equa- 
tion £Q|. As an example, if the Faraday effects are weaker than 
the absorption effect, i.e., n] 3> n v 2 and nj 3> Kg 2 , 
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is obtained to the leading order. The circular polarization is 
similarly given by ~ \rjv /t]i ~ k v/k;\. 

In the following sections, we apply this formulation to the 
late-time GRB afterglows. The anisotropic part of the di- 
electric tensor is tens of magnitude smaller than unity for the 
shocked plasma of a typical GRB afterglow. We assume that 
(1) the pitch-angle distribution of electrons is isotropic for 
simplicity; (2) the shocked plasma is spatially homogeneous; 
3 (3) the shocked plasma consists of a number of random cells 
within each of which magnetic field is ordered. With the third 
assumption, we obtain the observed linear and circular polar- 
izations by 1 / y/N times those for completely ordered mag- 
netic field, whe re N is the number of the random cells in the 
visibl e region dJones & O'Delll 119771 ; iGruzinov & Waxmanl 
1999). To reproduce the op tical detection at the level of 
~ 1 - 3% dCovino et alj|2004l) . N would be - 10 3 . 

3. POLARIZATION OF LATE-TIME GRB AFTERGLOWS 

In this section, we derive the polarization spectrum of the 
late-time afterglow, based on the standard external shock 
model in which all the electrons are accelerated, i.e., / = 1 
(see § 1). The energy distribution of the electrons is assumed 
to be dn/d-f e cx j~ p for j e > j m . The transfer coefficients 

3 Electron cooling makes the electron energy distribut ion inhomogene ous, 
but it can be neglected in the late phase of the afterglow I Sari et al. 1998). 




FIG. 1 . — Polarization spectrum of a late-time GRB afterglow, at t - 1 day, 
inferred from the standard external shock model in which all the electrons are 
accelerated, i.e., / = 1 . The degrees of linear polarization TIl (thin solid line) 
and circular polarization lie (thick solid line) are shown. The degrees are cal- 
culated as 1/ ~ 1CT 1 5 times those for completely ordered magne tic field, 
i.e., they are calibrated by detected optical linear polarization ( Co vino et all 
2004). Typical values of parameters are used: E = 10 s2 ergs, n = 1 cm -3 , 
e B = 10~ 2 , e e = 10-' , and p = 2.2. 

for such electr on plasma are summar ized for frequency re- 
gi on v > v m byjj ones AOIDelK 1977) and for 
by Matsumiva & Iokal (12003b . 4 where v m is the characteris- 
tic synchrotron frequency corresponding to j m and vb is the 
nonrelativistic electron Larmor frequency. 

The radius of the shock and the Lorentz factor of the 
shocked fluid evolve as R ~ (IJEt /AimipCri) 1 ^ and V ~ 
(17£' /10247rOT p c 5 »f 3 ) 1 / 8 , respectively, where t is the observer 
time dSariet alj|1998h . The comoving width of the shocked 
plasma shell can be estimated by R/4T, which we use as 
the path length of the transfer equation (Q~|). The mag- 
netic field strength, the minimum Lorentz factor and num- 
ber density of the accelerated electrons are written by B = 
(32inn p c 1 eBn) x / 1 T, j m = e e (m p /m e )T, and n acc = 4-Tn, respec- 
tively. Then we obtain v% : 
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Hz, respectively. Here (and 
hereafter) we have adopted the notation Q x = Q/\Q X in cgs 
units and td = t/\ day. 

Figure Q] illustrates the polarization spectrum of the late- 
time GRB afterglow. The frequencies at which r,7y, and tq 

equal unity are given by v a 
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Hz, where p = 2.2 has been used 
as a fiducial value. Since v a > vq ~ vy, so that no plasma ef- 
fects are significant in the optically thin regime v > v a and the 
intrinsic degree of polarization is obtained, LL: = i]q/t]i = 0.5 

and n c = \f]v/vi\ - 7m ("/"mY^ 3 for v < v m . For v > v m , 
n L = (P+1)/(P+ |) ^ 0.7 and n c ~ %t{v/v m y l l 2 . In the 
optically thick regime v < v a , t 2 3> Ty and t 2 ^> are satis- 
fied, the linear polarization is given by equation (f3J. Because 
Vq/Vi = k q/ K i = 0-5 for v <C v m , the intrinsic linear polariza- 
tion vanishes and n^, is only produced by the conversion of 
the circular polarization. The transfer equation dTJ indicates 
that I1 L w (K Q /n,)(r]v/ri, - Ky/m) «2x 10~ 2 \r}y /ru\ and 

^ ~ Tfy /rji — Ky / Kj ~ 6 x 10~ 2 \rjy /rjj\. All the characteristic 
frequencies v Bi vq, and vy are weakly dependent of time, so 
that the polarization spectrum does not evolve significantly. 

We adopt the value of = (a+ |)/2 different from that shown 



in Uone s & O'Dell (1977), and the sign of Ky should be changed in 
Matsumiva & Ioka 12003^ 
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The suppression of He due to absorption effect has not 
been pointed out in the context of GRB afterglows, since 
Matsumiva & Iokal d2003l) erroneously us ed Ky of opposite 
sign and lSagiv et al.1 (120041) neglected k v . ISagiv et al](l2004 
discussed similar propagation effects in early-time after- 
glows, but their discussion should be restricted to a frequency 
region v > v B (~ 10 10 - 10 11 Hz). For v > v B (> v a ), U L 
and lie they derived for the forward shock emissions are 
consistent with our results. 

4. A SIGNATURE OF THE THERMAL ELECTRONS 

Here we derive the polarization spectrum according to the 
standard external shock model in which the thermal elec- 
trons are left behind, i.e., / < 1 (see § 1), and show that 
the linear polarization may be suppressed even at frequen- 
cies higher than the absorption frequency v a . The elec- 
tron energy distribution is assumed to consist of the accel- 
erated electrons which are considered in § 3 and the ther- 
mal electrons with the Lorentz factor 7„, = T and the num- 
ber densit y n,h = 1(1- f)/f1>w, wh ere m e /m p < f < 1 (see 
Fig. 1 of lEichler & Wax man 2005). Then all the quanti- 
ties in this model can be written by using the parameters 
{E, n, eg, e e , p} as measured assuming / = 1, while the 
real values of the parameters are given by {£' = E/f, n 1 = 
n //> £ b = e #/> e e = e efi p' = p}- The characteristic syn- 
chrotron frequency of the thermal electrons is estimated by 

v m ~ 2 x lO^E^e^^C?'' 1 Hz. We approximate the trans- 
fer coefficients for the thermal electr ons as those for the mo- 
noenergetic distribution of electrons dSazonovll 1 969t iMelrosd 
Il980alfbh . Thus we consider the electron energy distribution 
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where K = (p- l)n acc 7 m 1 and H(x) is the Heaviside step func- 
tion. (Hereafter we describe the quantities related to the ther- 
mal electrons as Q.) For v ^S> v m , fjigy and K/,g.y damp expo- 
nentially. The remaining coefficients for the thermal electrons 
are different from those for the power-law distribution only by 
numerical factors. Here we show the expressions of Faraday 
coefficients, 
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where 9 is the angle between k and B and T^(x) is the Euler 
gamma function. The coefficients for the electron energy dis- 
tribution consisting of the thermal plus accelerated ones are 
given by the linear combination of the two contributions. 

In Figure [2] we show the linear polarization spectrum of 
the late-time GRB afterglow for the / < 1 model, compared 
with that for the / = 1 model obtained in § 3. The sensitivi- 
ties of ALMA and VLA for 1 hour integration time are also 
shown. They are derived by > F StU /F v , where F sl/ is the 
sensitivity limit for continuum radiation and F v is the flux of 
the afterglow. The flux is estimated by F v = F m (v/v m ) 1 ' 3 for 
v a < v < v m and F v = F m (i' a /i' m ) 1 ^(v/v a ) 2 for v m < v < v a , 

1-1 1/2 1/2 

where F m g 10 D 2 7 , 5^52^ £ b-2 m -ty anc ' ® * s t ' le luminosity 
distance (ISari et al.Ul998l) . 

For v ^> v m , the absorption effect of the thermal elec- 
trons is absent, and thus the absorption frequency is 
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FIG. 2. — Linear polarization spectra of a late-time GRB afterglow, at 
/ = 1 day, inferred from the standard external shock models in which all the 
electrons are accelerated (thin solid line; the same of 11/, in FigurefTJ and in 
which the thermal electrons are left behind with (1— f)/f = 10 (thick solid 
line). The parameters E, n, eg, e e , p., and N are the same as Figure[T] The 
dot-dashed lines describe the sensitivities of VLA and ALMA for a burst 
with D = 1 Gpc (z ~ 0.2) and an integration time of 1 hour. At frequencies 
lower than Oy, Hi is suppressed by the Faraday rotation effect of the thermal 
electrons, which can be detected with ALMA. 

the same as the / = 1 case. Since fy/ry — [(1 — 
/)//](7m/7m)" 2 (ln7m/ m 7m) > 1 ( see equation ©) and sim- 
ilarly tq/tq^S> 1, the Faraday effects are dominated by those 
of the thermal electrons. The ratio f v /fg ~ Qnj m )(v / v m ) 
is ^ 1 for small j m , i.e., at the late phase of the afterglow, so 
that the normal modes of this plasma is circularly polarized 
and the Faraday rotation effect is significant. The frequencies 
at which Ty and tq equal unity are given by 
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and %~4x 10 10 [(1- /)/ 10 fY'^E^nJ 3 e B 'l 2 t d 2li Hz, re- 
spectively. For v > vy, all the depths are smaller than unity, 
so that the intrinsic polarization is obtained. In the regime 
v a < v < vy, Ty ^S> tq ^> 1 ^> t is satisfied, so that LT^ is 
given by equation (|2j. It damps at low frequencies as cx v 2 
and oscillates with the period |Aj//j/| ~ 10~V 2 [. In the opti- 
cally thick regime v m <C v < v a , the transfer equation (JTJ for 
fy > tq > r > 1 indicat es that U L sa (fo/fy)(w / Vi~ K v/ K i) 
and lie ~ Vv /vi~ K v / ki d Jones & Q'Dell|[T977l) . and thus n L 
does not exceed He- For v i> m , both the absorption and 
the Faraday effects are dominated by the thermal electrons, 
and f 2 ^S> fy and f 2 ^ ¥q are satisfied. Then the polarization 
spectrum is similar to that for v < v a in the / = 1 model dis- 
cussed in § 3. It is important to note that both and lie are 
< 10~ 2 for v < v a , and they are far from detectable because 
the flux is suppressed in this regime (especially for v < v m , 
the additional absorption by the thermal electrons exists). 

The existence of the thermal electrons is characterized by 
the suppression of the linear polarization at v a < v < i>y. 
Necessary conditions for this suppression are vy ^> v a and 
vy ^S> vq. The former condition reduces to (1 — f)/f 3> 

Interestingly, the effect can be 
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seen even for as small number of thermal electrons as ( 1 - 
/) // ~ 1 _1 . The latter condition implies that t > 3 x 1 3 [( 1 - 

/) / ' f\i^ '£52^ '"o' 1 '' 25 ' e f/-2 s ' which shows that the late-time 
afterglow is suitable to search for the existence of the thermal 
electrons through the observation of linear polarization. 

If v a is determined by the observation of a bright burst and 
the linear polarization is not detected at v > v a with VLA 
and detected at v ^ v a with ALMA, it becomes clear that a 



4 



Toma et al. 



number of the thermal electrons exist and the magnetic field 
is ordered on large scales. If we determine vy, the electron- 
proton coupling parameter / can be constrained by equation 
0. 

5. DISCUSSION 

We have studied a signature of the thermal electrons only 
in the patchy coherent magnetic field model, while there 
are some other viable models for magnetic field configu- 
ration. In the mode l of random field with very short co - 
herence length (e.g., [Sarilll999t iGhisellini & Lazzati|[l999h . 
the coefficient k v averaged over the field configuration van- 
ishes, so that the Faraday de polarization of IT L does not 
occur dMatsumiva & Ioka 120031) . In the model of a com- 
binati on of random field B^ d and large-scale ordered field 
Z? or d dGranot & Koniglll2003l) . the depolarization by B or d can 
occur similarly as discussed in § 4. In this model ITl ~ 
0.7*^/(0 for v > v m and p = 2.2, so that B^/{B^ d ) ~ 
10 to reproduce the optical detection. Interestingly for 

Dy < v < v m n L = 0.52J*g/<fl*£) ~ 0.1. If such a high n L 
at vy < v <C v m is detected, it will be an evidence for the pres- 
ence of B or( j. 

Only upper limits have been obta ined so far for the radi o 
polarization from GRB afterglows dGranot & Taylor! 120051) . 
From brig ht GRB 030329, 11^ ~ 2% is measured in the op- 
tical band (iGreiner et al.l l2003). whereas 3<r limits < 1% are 
derived at 8.4 GHz. Such a low degree at radio may be at- 
tributed to the sourc e being optically t hick, since v a is esti- 
mated as ~ 19 GHz dTavlor et al.ll2005l) . 

If a large number of thermal electrons are left behind, i.e., 
/ < 1, the afterglow energy of GRBs should be E 1 = Ef~ l , 
where E is the afterglow energy estimated by / = 1 model 
and typical ly inferred to be ~ 10 51 5 ergs with jet collim ation 
correction (Pan aitescu & Kumdl20l)l lYost et al.ll2003l) . The 



association of GRBs with supernovae suggests that / > 10~ L5 
is a conservative lower limit. The energy of prompt 7-ray 
emission is typically similar to E, and some of the models 
of early-time afterglows i mply that the efficiency of the 7-ray 
emission is > 90% (e.g., Iloka et al]l2006t iToma et al.ll2006l: 
iGranot et alll2006t iFan & Pirar]l2006l) . If the external shock 
model with / < 1 is applicable to early-time afterglows, the 
7 -ray efficiency proble m would be solved. 

iMundell et al.1 (120071) have reported a 2a upper limit < 8% 
on the optical polarization in early-time afterglow of GRB 
060418 (t ~ 200 s), and argued that presence of a large-scale 
ordered field in the GRB jet is ruled out. However, the ro- 
tation frequency vy for a typical reverse shocked ejecta with 
ordered field is ~ 10 15 Hz, so that the low level of polariza- 
tion degree would result from the Faraday depolarization (see 
ISagiv et al.f2004l) . 

It is suggested that there are electrons well coupled to 
protons at som e nonrelativistic collisionless shocks, that is, 
e e ~ 0.5 (e.g.. iMarkevitchlfeOOfSl) . However, the fraction / of 
total electrons that are coup led to protons has not been dis- 
cussed seriously. Recently, ISpitkovskvl d2007l) has reported 
the results of 2-dimensional particle-in-cell simulations of rel- 
ativistic collisionless shocks in electron-proton plasma with 
realistic value of m p /m e that e e w 0.5 and / ~ 1 were re- 
alized. However, it is too early to interpret the results con- 
clusive, since long-term 3-dimensional simulations with good 
resolution have not been done and it has not been understood 
whether the heated electrons are accelerated into the power- 
law energy spectrum. The late-time radio polarimetry may be 
an important test for more realistic simulations and theories. 
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